Measurement of D and meson masses with the KEDR Detector 
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, Abstract 

■The masses of the neutral and charged D mesons have been measured with the KEDR detector at the VEPP-4M electron-positron 
coUider: 

Mflo =1865.30 + 0.33 + 0.23 MeV, 
Md- =1869.53 + 0.49 + 0.20 MeV. 

'Key words: D meson, charm, mass, X(3872), ifi{377Q) 
•PACS: 13.20.Fc, 13.20.Gd, 14.40.Lb 
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1. Introduction 

Neutral and charged D mesons are the ground states in the 
'family of open charm mesons. Measurement of their masses 
provides a mass scale for the heavier excited states. In addition, 
a precise measurement of the meson mass should help to 
understand the nature of the narrow X(3872) state idHSlil, 
which, according to some models, is a bound state of £)" and 
(D*^ mesons Jst] and has a mass very close to the sum of the 
D° and D'° meson masses. Presently, the world-average D 
mass value fSI (M^o = 1864.84 + 0.17 MeV) is dominated 
by the CLEG measurement Moo = 1864.847 + 0.150(stat) + 
0.095(syst) MeV fl, which uses the decay D° 4>K'^. Other 
D meson mass measurements are much less precise. These 
measurements were carried out long ago in the MARK-II ex- 
periment at the SPEAR e+e" collider Istl, and by the ACC- 
MOR collaboration in a fixed-target experiment jsll. Both mea- 
surements are dominated by the systematic uncertainty, which 
in the case of MARK-II is related to beam energy calibra- 
tion. In addition, the mass of the is constrained by the 
D" mass and a mass difference Mo* - M^o much more pre- 
cisely than directly measured: the world-average mass is 
Md+ = 1869.62 + 0.20 MeV, while the direct measurements 



yield Mo+ = 1869.5 + 0.5 MeV. 

As both D° and mass values are based on a single mea- 
surement, the cross-check involving a method different from the 
one used at CLEO is essential. This paper describes a measure- 
ment which has been performed with the KEDR detector at the 
VEPP-4M e^e- collider using the decay (/r(3770) DD. 



2. Experimental facility 

The electron-positron accelerator complex VEPP-4M ifiol] 
designed for high-energy physics experiments in the center-of- 
mass (CM) energy range from 2 to 12 GeV is currently running 
in the iff family region. The collider consists of two half -rings, 
an experimental section where the KEDR detector is installed, 
and a straight section, which includes an RF cavity and injec- 
tion system. The circumference of the VEPP-4M ring is 366 m. 
The luminosity at the J/i// in an operation mode with 2 by 2 
bunches reaches X. = 10^*' cm"-s"'. 

Precise measurement of beam energy can be performed 
at VEPP-4M using the resonant depolarization method ll 111 . 
The method is based on the measurement of the spin preces- 
sion frequency of the polarized beam, which depends on its 
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energy. Using resonant depolarization, the precision of the 
beam energy rneasurement reached in the KEDR experiment 
is^lOkeVlH. 



The KEDR detector 11 1311 includes a tracking system consist- 
ing of a vertex detector and a drift chamber, a particle identifica- 
tion (PID) system of aerogel Cherenkov counters and scintilla- 
tion time-of-flight counters, and an electromagnetic calorime- 
ter based on liquid krypton (in the barrel part) and Csl crys- 
tals (endcap part). The superconducting solenoid provides a 
longitudinal magnetic field of 0.6 T. A muon system is in- 
stalled inside the magnet yoke. The detector also includes a 
high-resolution tagging system for studies of two-photon pro- 
cesses. The online luminosity measurement is performed with 
sampling calorimeters which detect photons from the process 
of single brehmsstrahlung. 

Charged tracks are reconstructed in the drift chamber (DC) 
and vertex detector (VD). DC 11411 has a cylindrical shape of 
1 100 mm length, an outer radius of 535 mm and is filled with 
pure dimethyl ether DC cells form seven concentric layers: 
four axial layers and three stereo-layers to measure track coor- 
dinates along the beam axis. The coordinate resolution aver- 
aged over drift length is 100 jum. VD llisll is installed between 
the vacuum chamber and DC and increases a solid angle ac- 
cessible to the tracking system to 98%. VD consists of 312 
cylindrical drift tubes aligned in 6 layers. It is filled with an 
Ar-l-30%CO2 gas mixture and has a coordinate resolution of 
250 jum. The momentum resolution of the tracking system is 
o-pip = 2% e (4% X p[GeV]). 

Scintillation counters of the time-of-flight system (TOP) are 
used in a fast charged trigger and for identification of the 
charged particles by their flight time. The TOP system con- 
sists of 32 plastic scintillation counters in the barrel part and in 
each of the endcaps. The flight time resolution is about 350 ps, 
which corresponds to njK separation at the level of more than 
two standard deviations for momenta up to 650 MeV. 

Aerogel Cherenkov counters (ACC) [Ifijo are used for particle 
identification in the momentum region not covered by the TOP 
system and ionizations measurements in DC. ACC uses aero- 
gel with the refractive index of 1.05 and wavelength shifters for 
light collection. This allows one to identify n and K mesons 
in the momentum range of 0.6 to 1.5 GeV The system design 
includes 160 counters in the endcap and barrel parts, each ar- 
ranged in two layers. During data taking only one layer of ACC 
was installed, and it was not used because of insufficient effi- 
ciency. 

The barrel part of the electromagnetic calorimeter is a liquid 
krypton ionization detector UtIi . The calorimeter provides an 
energy resolution of 3.0% at the energy of 1.8 GeV and a spa- 
tial resolution of 0.6-1 .0 mm for charged particles and photons. 
The endcap part of the calorimeter is based on 1536 CsI(Na) 
scintillation crystals jfisf] with an energy resolution of 3.5% at 
1.8 GeV, and a spatial resolution of 8 mm. 

The muon system lfl9ll is used to identify muons by their 
flight path in the dense medium of the magnetic yoke. It con- 
sists of three layers of streamer tubes with 74% solid angle cov- 
erage, the total number of channels is 544. The average longi- 
tudinal resolution is 3.5 cm, and the detection efficiency for the 



most of the covered angles is 99%. 

Trigger of the KEDR detector consists of two levels: pri- 
mary (PT) and secondary (ST). Both PT and ST operate at the 
hardware level. PT uses signals from TOP counters and both 
calorimeters as inputs, the typical rate is 5 10 kHz. ST uses 
signals from VD, DC and muon system in addition to systems 
listed above, and the rate is 50 H- 150 Hz. 



3. Measurement method 

Measurement of D meson masses is performed using the 
near-threshold e^e^ — > DD production with full reconstruction 
of one of the D mesons. Neutral D mesons are reconstructed 
in the K^n^ final state, charged D mesons are reconstructed 
in the K^n^n^ final state (charge-conjugate states are implied 
throughout this paper). To increase a data sample, the collider is 
operated at the peak of the i/r(3770) resonance. The production 
cross sections at this energy are cr(D"D") - 3. 66+0. 03+0. 06 nb 
and o-(D^D-) = 2.91 + 0.03 + 0.05 nb (Q. 

The invariant mass of the D meson can be calculated as 
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(so-called beam-constrained mass), where fibeam is the average 
energy of colliding beams, pi are the momenta of the D decay 
products. The mass calculated this way is determined more pre- 
cisely than in the case when the D energy is obtained from the 
energies of the decay products. The precision of Mo measure- 
ment in one event is 
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where c^^' is the CM energy spread. The contribution of the 
momentum resolution is suppressed significantly due to small 
D momentum (po - 260 MeV). 

In addition to Mbc, D mesons are effectively selected by the 
CM energy difference 
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Mf +pf- £be: 
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where M,- and p, are the masses and momenta of the D decay 
products. The signal events should satisfy a condition AE ^ 
0. In our analysis, we select a relatively wide region of Mbc 
and AE close to Mbc ~ Mo and AE ~ (specifically, Mbc > 
1700 MeV, lA^I < 300 MeV); then a fit of the event density 
is performed with D mass as one of the parameters, with the 
background contribution taken into account. The background 
in our analysis comes from the random combinations of tracks 
of the continuum process e^e^ qq {q - u, d, s), from other 
decays of D mesons, and from the signal decays where some 
tracks are picked up from the decay of the other D meson. 

While calculating Mbc, we employ a kinematic fit with the 
AE - constraint. It is done by minimizing the function 
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formed by the momenta of the daughter particles 

ip'i-Pif 



X 



= 2 



(4) 



where p, and cTp, are the measured momenta of the daugh- 
ter particles and their errors obtained from the track fit, re- 
spectively, and p'. are the fitted momenta which satisfy the 
AE{p'j) - constraint. The use of Mbc constructed from the fit- 
ted momenta results in a certain improvement of its resolution 
and significantly reduces the dependence of measured mass on 
the absolute momentum calibration (see below). 

The precision of the momentum measurement has direct in- 
fluence on the D mass measurement. The following sources of 
momentum reconstruction uncertainties are considered in our 
analysis: 

1 . Simulation of ionization losses in the detector material. 

Reconstruction of cosmic tracks is used to check the va- 
lidity of the simulation. We select the cosmic tracks that 
traverse the vacuum chamber and fit their upper and lower 
parts separately. The average difference of the upper and 
the lower track momenta due to energy loss in the detector 
material is compared with the result of the simulation. 

2. Absolute momentum calibration (this is equivalent to 
the knowledge of the average magnetic field in the tracking 
system), described by the scale coefficient a which relates 
the true track momentum paue and the measured momen- 
tum p: 

Ptiue = ap . (5) 

Then 
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The momentum scale can be calibrated using the same 
events as in the D mass measurement by measuring the 
average bias of the AE value: 
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Sensitivities to the scale coefficient a are given by 

dAE pi pi 

^ + i2L ^ 1580 MeV 

da Ek E„ 



for D" K decay, and 
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da 



^ 1490 MeV 



(8) 
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for K^ji^n^ decay. The numerical values of the 

dAE I da derivatives are obtained using the Monte Carlo 
(MC) simulation of the corresponding decays. When 
the kinematic fit with AE = is employed for an Mbc 
calculation, such a correction is effectively applied to 



each event, and thus the dependence of Mo on the ab- 
solute momentum calibration is significantly reduced (to 
dMo/da = -3 MeV for D" -> K-7T+ and -12 MeV for 
D+ ^ K-n^7T+). 

As a cross-check, we also use other processes for the ab- 
solute momentum calibration: the inclusive A'^' — > tt'^tt^ 
reconstruction and e^e^ 'I'i^S) J/t/zn^n^ process. 
3. Simulation of the momentum resolution. Since the D 
meson sample is limited, we use full MC simulation of 
the detector to determine the shapes of the signal distribu- 
tions. The description of the momentum resolution in the 
simulation is adjusted using events of elastic e^e^ scatter- 
ing, inclusive resonstruction of K'^ — > n'^n^ decay, and the 
process e"^e" (^(25) — > Jltj/n^Ti^. 

4. Analysis of /)" -» K~n^ 

The analysis uses a sample of 0.9 pb ' accumulated with the 
KEDR detector at the energy of the i/'(3770) resonance. Mul- 
tihadron candidates which contain at least three tracks close 
to the interaction region (transverse distance from the beam 
R < 5 mm, and longitudinal distance \z\ < 120 mm) forming 
a common vertex are selected at the first stage of the analysis. 
The pairs of oppositely charged tracks are taken as decay 
candidates with the following requirements: 

• Number of track hits A^hits ^ 24 (the maximum number of 
hits per track is 48), 

• Track fit quality < 50, 

• Transverse momentum: 100 MeV < pj < 2000 MeV. 

• Energy of the associated cluster in the calorimeter E < 
1000 MeV. 

We expect around 100 K^n^ signal events for this 

sample. In order to measure the mass most efficiently, the 
unbinned maximum likelihood fit procedure is used. Except 
for the Mbc variable, the likelihood function includes two other 
variables which allow one to efficiently separate the signal from 
the background: the energy difference AE - Ed - £beam ® and 
the difference of the absolute values of momenta for D decay 
products in the CM frame A\p\. 

The UkeUhood function has the form: 

N „ 

- 2 log Z{a) = -2 ^ log p(Vi|a) + 2N log I p(\\a)d\ ,(11) 

!=0 ^ 

where v = (Mbc, AE, A|/?|) are the variables that character- 
ize one event, pi\\a) is the probability distribution function 
(PDF) of these variables depending on the fit parameters a - 

iMD,{AE),b,d„bDDy- 

p(\\a) = Psig(y\MD, (AE)) + b„dsPuds{y) + boDPooiy) ■ (12) 

Here p^tg is the PDF of the signal events which depends on 
Mo (D" mass) and (AE) (the central value of the AE distribu- 
tion), p„ds is the PDF for the background process e^e^ qq 
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(q - u, d, s), and pon is the PDF for the background from 
e^e" — > DD decays with D decaying to all modes other than 
the signal one, buds and boo are their relative magnitudes. The 
shape of the psig, p„ds and poo distributions is obtained from 
the MC simulation. Such a fit procedure gives only a shape of 
the fitted distribution without the absolute normalization. The 
numbers of signal and background events can be extracted by 
taking the total number of events in a sample and fractions of 
the corresponding events from the fit. 

For a proper calculation of AE - E„ + Ek - £beam, the k/K 
identification is needed. Presently it cannot be performed reli- 
ably in the momentum range near 800 MeV. Fortunately, since 
the D meson momentum is small, the momenta of K and n dif- 
fer by a small amount, and the maxumum error (in the case of 
wrong mass assignment) is not larger than 30 MeV. Thus, we 
take the following combination as a D meson energy: 



where 



E' = (Ek-„* + Ek*„-)/2 , 



(13) 



(14) 



The energy E' calculated this way is practically unbiased from 
the true energy E. A bias can appear if the detection efficiency 
varies with momenta of the final state particles; we estimate the 
upper limit of this bias to be 1 .5 MeV. According to ^ and (|9]i, 
this bias is propagated to an M^o bias of 0.034 MeV. E' differs 
from E by less than 15 MeV in each event, this only slightly 
affects the AE error due to momentum resolution. 

Use of the /S\p\ variable allows us to obtain an estimate of the 
Mbc resolution on the event-by-event basis, thus improving the 
overall statistical accuracy of the measurement. We use the fact 
that this resolution depends strongly on decay kinematics — it 
can be up to three times better for events where the daughter 
particles from decay move transversely to the direction of 
the Z)" (Alpl is around zero for these events), than for events 
where they move along this direction (see Fig.[T]l- 

The variables Mbc and A|p| use the momenta of the daugh- 
ter particles after the kinematic fit with the AE - constraint, 
while AE is calculated using uncorrected momenta. We select 
combinations that satisfy the following requirements for the fur- 
ther analysis: Mbc > 1700 MeV, lA^I < 300 MeV. 

Simulation of signal events is performed with the MC gener- 
ator for e^e" — > DD decays where D-meson decays are simu- 
lated by the JETSET 7.4 package and the radiative cor- 
rections are taken into account in both initial (ISR, using the 
RADCOR package IXM with Kuraev-Fadin model I23il ). and final 
states (FSR, the PHOTOS package Q)- The ISR corrections 
use the e^e^ — > DD cross section dependence of the resonant 
production of the iff{3nQ) according to a Breit-Wigner ampli- 
tude with M = 3771 MeV and F = 23 MeV Q, without the 
nonresonant contribution and taking into account phase space 
dependence at the production threshold. The full simulation of 
the KEDR detector is performed using the GEANT 3.21 pack- 
age HI]. 
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Figure 1: Correlation of M\,^ and A\p\ variables for ^ K n"^ decays (MC 
simulation). 



The PDF of the signal events psig is a function of three pa- 
rameters Mbc, AE, and A|p|. It is parameterized with the sum 
of two two-dimensional Gaussian distributions in Mbc and AE 
(representing the core and the tails of the distribution) with a 
correlation and with the quadratic dependence of the Mbc res- 
olution on A\p\. The core distribution is asymmetric in Mbc 
(with the resolutions cri(Mbc) and crR(Mbc) for the left and right 
slopes, respectively). The A\p\ distribution is uniform with a 
small quadratic correction and with the kinematic constraint 
(A|/7|)^ < - M^. The parameters of the signal distri- 

bution are obtained from the fit to the simulated signal sam- 
ple. The core resolutions obtained from the MC for Mbc are 
o-z.(Mbc) = 0.98 + 0.03 MeV, crR(Mbc) = 2.45 + 0.06 MeV 
(at A\p\ - 0), the Mbc resolution at A\p\ - 200 MeV is 
4.6 + 0.1 MeV, the core resolution of AE is 48.3 + 0.3 MeV. 

The background from the continuum e^e" qq process 
(where q - u, d, s) is simulated using the JETSET 7 . 4 e^e" — > 
qq generator. The PDF is parameterized as 



PudAMhc, AE, A\p\) = exp 



El 

beam 



■ k2AE 



(15) 



{\+hA\p\^), 



where ki are free parameters. The kinematic limit at Mbc = 
-Ebeam IS provided by the (Alpl)^ < EI^^^ - Mg^ constraint. 

The background from e^e" — » DD decays is simulated us- 
ing the JETSET 7 . 4 generator, where the signal process Z)" — > 
K^Ti^ is suppressed in the decay table. The PDF for DD back- 
ground is parameterized with the function poo of the same form 
as for Puds, with the addition of three two-dimensional Gaus- 
sian distributions in Mbc and AE. Two of them describe the 
background from D^ — > n^n'' and D" — > K^K , while the third 
one is responsible for the decays of D mesons to three and more 
particles. 

The combinatorial background coming from the signal 
events where one or more tracks were taken from the decay 
of the other D meson, were studied using the signal MC sam- 
ple. The distribution of fit variables for these events is similar to 
the background from the continuum events, and their fraction is 
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Table 1: Results of the fit to the D° ^ R-n* data sample Table 2: Results of the fit to the D* ^ K tt+k* data sample 



Md 


1865.05 ± 0.33 MeV 




Mo 


1869.58 + 0.49 MeV 


(AE) 


-0.7 + 7.3 MeV 




(AE) 


2.5 + 5.0 MeV 


Number of signal events 


98.4+13.1 




Number of signal events 


109.8 + 15.3 


Number of qq events 


18.3 + 2.4 




Number of qq events 


85.3 + 11.8 


Number of DD events 


4.8 + 0.8 




Number of DD events 


11.4 + 2.2 



2.5% of the number of signal events, which is negligible com- 
pared to the continuum contribution. We therefore do not treat 
this background separately, and its contribution is effectively 
taken by the continuum component. 

The result of the fit to the experimental data is shown in 
Fig. |2l In the fit we use the function ( fT2] ) with Md, (AE) as 
well as the relative magnitudes of the continuum and DD back- 
grounds as free parameters. 

The momentum correction coefficient a is chosen to keep the 
value of {AE} close to zero. Event selection is performed with 
a - 1 .030; after the residual AE bias is taken into account its 
value is ff = 1.0304 ± 0.0046. The results of the fit are shown 
in Table [T] The numbers of events are presented for the signal 
region |A£| < 100 MeV, 1855 MeV< Mbc < 1875 MeV. 

To obtain the D" mass, one has to take into account a possible 
deviation of the fit parameters Mo and (AE) from the true D^ 
mass and energy. In particular, the central value of Mo can be 
shifted due to the asymmetric resolution function and radiative 
corrections. This deviation is corrected using the MC simu- 
lation. The final value of the D° mass after the correction is 
Moo = 1865.30 + 0.33 MeV 

5. Analysis of -» K'n'^n'* 

The three-body decay D^ — > K'n^n^ has more kinematic 
parameters and there is no simple variable (such as Alpl in the 
D" — » K n^ case), which determines the precision of the Mbc 
reconstruction. Therefore, we use only two variables, Mbc and 
AE, in a fit of this mode. 

The mode D^ —> K^n^n^ does not have a problem with njK 
identification for the AE calculation, since the sign of the kaon 
charge is opposite to the pion charges and thus energies of all 
the particles can be obtained unambiguously. The triplets of 
tracks with the charge of one of the tracks ("kaon") opposite to 
the charges of the two other tracks ("pions") are taken as D"^ 
decay candidates. 

The requirements for the track selection are the same as in 
the D" — > K n^ case. Since the significant part of the kaon 
tracks in the three-body decay have relatively low momen- 
tum (under 500 MeV), an additional suppression of the back- 
ground from pions is possible using the TOP system. The se- 
lection uses the following requirement on the flight time for a 
kaon candidate, which hits the barrel part of the TOP system: 
ATjoF - TjoF - Tk(pk) > -0.8 ns (or 2.3 times the flight 
time resolution), where Tk{pk) is the expected flight time for a 
kaon with the momentum px and Ttof is the measured flight 
time. As a result of this requirement the background fraction 



is reduced by a factor of 2.3 for the continuum background and 
3.3 for the DD background. 

The Mbc variable uses the momenta of the daughter particles 
after the kinematic fit with the AE - constraint. The variable 
AE is calculated using uncorrected momenta. We select com- 
binations that satisfy the following requirements for the further 
analysis: Mbc > 1700 MeV, |A£| < 300 MeV. 

As in the case of D^ — > K^n^ decay, simulation is performed 
using the e^e^ DD generator taking into account the ISR 
and PSR effects. The signal PDP psig is parameterized in the 
same way as for the D^ K^n^ mode, but without A|p| depen- 
dence. The core resolutions obtained from the MC for Mbc are 
a-i.(Mbc) = 2.07 + 0.05 MeV, cr«(Mbc) = 2.52 + 0.06 MeV, the 
core resolution of AE is 26.5 + 0.4 MeV. 

To parameterize the continuum e^e^ — > qq background, we 
use the empirical function of Mbc proposed in the Argus exper- 
iment 12711 and the exponent of the quadratic form in AE: 



PudAMbc, AE) = yexp [-kiy' - [ka + k^y^jAE + k^AE^) , 

(16) 

where y — -^1 - A/^c/^beam- '^^^ coefficients ki are free pa- 
rameters in the fit. The coefficient k^, is responsible for the Mbc 
dependence of the AE slope, which appears after the kinematic 
fit to AE = 0. The PDP for the e+e^ -> DD background poo 
is parameterized with the distribution of the same form as for 
Puds, with the addition of two two-dimensional Gaussian dis- 
tributions in Mbc and AE. They describe the contributions of 
D^ K^K^n^ and D decays to four and more particles. The 
combinatorial background from the signal events as in the case 
of the D" — > K^ji^ mode is effectively taken into account by 
the continuum component. 

The result of the fit to the data is shown in Pig. [3] The mo- 
mentum correction factor a is chosen such that {AE) is close 
to zero. The value a - 1.027 is used for event selection, 
and after taking into account the residual AE bias its value is 
a - 1.0252 + 0.0035. The results of the fit are shown in Ta- 
ble 12] The numbers of events are shown for the signal region 
|A£| < 70 MeV, 1860 MeV< Mbc < 1880 MeV. 

As in the case of the D" — > K^ji^ mode, the D^ mass ob- 
tained in the fit is corrected for the bias of Mo and AE using 
MC simulation. The value of the D^ mass after the correction 
is Mo* = 1869.53 + 0.49 MeV. 

6. Systematic uncertainties 

The estimates of systematic uncertainties in the D mass mea- 
surements are shown in Table [3] 
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Figure 2: Experimental data (points with the error bars) and the results of the fit (histogram) for the ^ K 71^ decay. Mbc distribution for events with 
|A£:| < 100 MeV (a), AE distribution for events with 1855 MeV< < 1875 MeV (b), and the experimental (Mt,^,AE) scatter plot (c). 
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Figure 3: Experimental data (points with the error bars) and the results of the fit (histogram) for the D* — > K n'^'n'*' decay. Mbc distribution for events with 
|A£| < 70 MeV (a), AE distribution for events with 1860 MeV< < 1880 MeV (b), and the experimental (Mbc, A£) scatter plot (c). 



Table 3: Systematic uncertainties in the and mass measurements 





AMflO, MeV 


AMfl+, MeV 


Absolute momentum calibration 


0.04 


0.04 


Ionization loss in material 


0.01 


0.03 


Momentum resolution 


0.13 


0.10 


ISR corrections 


0.16 


0.11 


Signal PDF 


0.07 


0.05 


Continuum background PDF 


0.04 


0.09 


DD background PDF 


0.03 


0.06 


Beam energy calibration 


0.01 


0.01 


Total 


0.23 


0.20 



The contribution of absolute momentum calibration is deter- 
mined by the precision of the {AE} measurement and is prop- 
agated to the uncertainty of the mass measurement using the 
c/Mbc/iia dependence. For the £)" — > K^tt^ mode, the addi- 
tional factor, which dominates the momentum calibration error, 
is a possible bias of the approximate D energy calculation using 
Eq. ( fT3l l in the absence of tt/K identification. However, due to 
smaller dM\,c/da value the momentum calibration uncertainty 
for this mode is close to the one for the K^n^n^ mode. 

The uncertainty of the simulation of ionization losses in the 
detector material is estimated by the variation of the corre- 
sponding correction term within the limits given by the cosmic 
track measurement (+20%). 

The uncertainty due to momentum resolution is estimated by 
using different procedures matching the resolution in the simu- 
lation with the experimental one (either by introducing a cor- 
rection to the drift curve of the DC, or by smearing the re- 
constructed momenta) and by varying the tuning parameters 
responsible for the momentum resolution matching within the 
limits given by the calibration processes. 

The ISR correction uncertainty is dominated by the uncer- 
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tainty of the energy dependence of the cross section cr(e+e — > 
DD). The defauh fit uses the iA(3770) parameters from PDG- 
2006 for the cross section (M = 3771.1 + 2.4 MeV and T = 
23.0 + 2.7 MeV 12511 ). To estimate a systematic error, these 
parameters are varied within their errors, also the PDG-2008 
value is used (F - 27 MeV |@]). In addition, the nonresonant 
contribution is added incoherently to the 1 nb cross section at 
the peak. The quadratic sum of deviations in M, T 

and non-resonant contribution is taken as the systematic error. 
The model uncertainty of Kuraev-Fadin formulae UM is small 
(~ 0.1%) and has a negligible effect on our results. 

The uncertainty due to signal shape parameterization is esti- 
mated by using the alternative shape with one Gaussian peak. 

The continuum background shape uncertainty is estimated 
by using the alternative generator for the system of pions with 
the varying multiplicity in the simulation, and also by relax- 
ing the background shape parameters in the experimental fit. 
The contribution of the DD background shape is estimated by 
relaxing the relative magnitude of the Gaussian peaks and the 
non-peaking component in the experimental fit, and by exclud- 
ing one of the Gaussian peaks from the background shape pa- 
rameterization. In the case of the — > K'n^n^ mode, the 
background shape variation also includes the shapes obtained 
without a TOF requirement to take into account the uncertainty 
in the TOF simulation. 

To check possible inconsistencies in the three-dimensional 
signal and background description of the D" — > K n^ mode, we 
perform separate fits to data with different A\p\ requirements. 
The results are consistent within statistical errors. 

The error of the beam energy calibration is dominated by the 
precision of the beam energy interpolation between successive 
energy measurements using the resonant depolarization tech- 
nique. It does not exceed 70 keV and is of order 10 keV for 
most of the data sample. The uncertainty due to beam energy 
calibration is estimated in the worst case of a 100% correlation 
between all energy measurements. 



7. Conclusion 



MARK II 



ACCMOR 



CLEO 



KEDR 



PDG: 1864.84± 0.17 
tO.60 



186' . 85 ± 0.21 



1865.30+0.40 



1864 1865 

D° mass (MeV) 



MARK II 



ACCMOR 



KEDR 



1869.40±0.60 



PDG: 1869.62 ± 0.20 



1 870.00 ± 1.12 



186!.53±0.53 



1869 1870 1871 

D"^ mass (MeV) 



Figure 4: Compaiison of D meson masses with the other measurements. The 
thiclc and thin eiTor bai's show the systematic and the total errors, respectively. 
The shaded areas are the PDG-2008 values @|. The PDG value for the is ob- 
tained using the measured mass difference of the and Z)" mesons. MARK-II 
does not quote the systematic error separately. 



Masses of the neutral and charged D mesons have been mea- 
sured with the KEDR detector at the VEPP-4M e+e" collider 
operated in the region of the 1/^(3770) meson. The analysis uses 
a data sample of 0.9 pb ' with D mesons reconstructed in the 
decays D° — > K n^ and D^ — > K'n^n^ . The values of the 
masses obtained are 

• Moo = 1865.30 + 0.33 + 0.23 MeV, 

• Md- = 1869.53 + 0.49 + 0.20 MeV. 

The D° mass value is consistent with the more precise measure- 
ment of the CLEO collaboration 101, while that of the D^ mass 
is presently the most precise direct determination. 

Comparison of the D meson masses obtained in this analysis 
with the other measurements is shown in Fig.|4l 
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